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The statistical mechanical description of two-dimensional inviscid fluid turbulence is re- 
considered. Using this description, we make predictions about turbulent flow in a rapidly 
rotating laboratory annulus. Measurements on the continuously forced, weakly dissipative 
flow reveal coherent vortices in a mean zonal flow. Statistical mechanics has two crucial 
requirements for equilibrium: statistical independence of macro-cells (subsystems) and 
additivity of invariants of macro-cells. We investigate these requirements in the context 
of the annulus experiment. The energy invariant, an extensive quantity, should thus be 
additive; i.e., the interaction energy between a macro-cell and the rest of the system 
(reservoir) should be small, and this is verified experimentally. Similarly, we use ad- 
ditivity to select the appropriate Casimir invariants from the infinite set available in 
vortex dynamics, and we do this in such a way that the exchange of micro-cells within 
a macro-cell does not alter an invariant of a macro-cell. A novel feature of the present 
study is our choice of macro-cells, which are continuous phase space curves based on 
mean values of the streamfunction. Quantities such as energy and enstrophy can be de- 
fined on each curve, and these lead to a local canonical distribution that is also defined 
on each curve. The distribution obtained describes the anisotropic and inhomogeneous 
properties of a flow. Our approach leads to the prediction that on a mean streamfunction 
curve there should be a linear relation between the ensemble-averaged potential vorticity 
and the time-averaged streamfunction, and our laboratory data are in good accord with 
this prediction. Further, the approach predicts that although the probability distribu- 
tion function for potential vorticity in the entire system is non-Gaussian, the distribution 
function of micro-cells should be Gaussian on the macro-cells, i.e., for curves defined by 
mean values of the streamfunction. This prediction is also supported by the data. While 
the statistical mechanics approach used was motivated by and applied to experiments 
on turbulence in a rotating annulus, the approach is quite general and is applicable to a 
large class of Hamiltonian systems, including drift-wave plasma models, Vlasov-Poisson 
dynamics, and kinetic theories of stellar dynamics. 



1. Introduction 

1.1. Overview 

Statistical mechanics provides a way to calculate the macroscopic properties of matter 
from the behavior of microscopic constituents. Instead of considering all motions of the 
individual constituents, one describes observable quantities averaged over constituent 
Hamiltonian trajectories, and averages are evaluated using the probability distribution 
of possible microstates. Likewise, fluid systems with a local balance between dissipation 
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and forcing have been described by statistical mechanics with the inclusion of constraints 
based on invariants of the dynamics. In general, such statistical theories for fluids are 
based on the idea that the macroscopic behavior of the fluid turbulence can be described 
without knowing detailed information about small scale vortices. 

The justification of statistical mechanics based on ideal two-dimensional fluid equations 
is open to question, given the existence of forcing, dissipation, three-dimensional effects, 
temperature gradients, etc., that certainly occur in real fluid flows. Moreover, one must 
square the idea of cascading with the approach to statistical equilibrium. Ultimately such 
a justification is very difficult and would rely on delicate mathematical limits. However, 
its success amounts to the idea that the fluid system can in some sense be described by 
weakly interacting subsystems, where the behavior of a single subsystem can be described 
by weak coupling to a heat bath that embodies all of the other subsystems and all of the 
omitted effects. In the end 'the proof of the pudding is in the eating', and our justification 
is based on experimental observations. 

Intimately related to the existence of subsystems is the question of which invariants 
to incorporate into a statistical mechanics treatment of fluids. One aim of the present 
paper is to investigate this question. We investigate this question both conceptually 
and experimentally and come to the conclusion that quadratic invariants (energy and 
enstrophy) are most important. Our conclusion follows from the observation that these 
invariants possess the property of additivity. 

The microscopic dynamics of conventional statistical mechanics is finite dimensional, 
but to describe macrosopic phenomena one takes the thermodynamic limit in which 
the number of degrees of freedom tends to infinity. However, the dynamics of a two- 
dimensional fluid is already infinite dimensional and possesses an infinite number of 
invariants; so, in order to make progress with a statistical mechanics approach one must 
extract a finite-dimensional model, and such a model cannot conserve all of the in- 
variants of the original fluid system. In calculations one may also take limits of this 
finite-dimensional model, but the results of these limits may depend upon which of the 
invariants are maintained. Additivity of macroscopic invariants and statistical indepen- 
dence of subsystems are crucial properties in conventional statistical mechanics [see e.g. 
Landau & Lifshitz (1980)]. Because not all invariants of a system are additive, this prop- 
erty can be used to select invariants for statistical mechanics from the infinite number 
possessed by two-dimensional fluid systems. 

Related to the choice of additive invariants is the choice of subsystems. This choice 
requires the identification of two scales, a macroscopic scale and a microscopic scale, 
which we call A and 5, respectively, and phase space cells of these characteristic sizes are 
considered. In classical statistical mechanics, the micro-cells usually refer to individual 
particles, while the macro-cells, the subsystems, are selected to be large enough to contain 
many particles yet small enough to have uniform invariants. We address in detail the 
choice of these cells for the fluid in §5, but it is clear that a macro-cell should contain many 
micro-cells, yet be small enough so that the vorticity and streamfunction are constant. 
This condition is sufficient for statistical independence, but the converse is not always 
true. In any event we seek to define macro-cells that are nearly statistically independent 
and consider only invariants that are additive over these cells. 

A second aim of the present work is to propose the idea that temporal mean values 
of the streamfunction provide a natural coordinate system for describing inhomogeneous 
turbulence, a coordinate system that can be used to define statistically independent sub- 
systems. We suggest this idea because contours of the streamfunction for two-dimensional 
inviscid fluid flow tend to be smooth and because there tends to be a strong statistical 
dependence of vorticity or potential vorticity along those contours. Streamfunction con- 
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tours are much smoother than vorticity contours because of the smoothing property of 
the inverse Laplacian. Therefore, there is a natural separation of length scales: the large 
scale associated with variation of the streamfunction contours and the fine scale that is 
needed to resolve the vorticity or potential vorticity. We take these to be our scales A 
and 5, respectively. We test this idea experimentally by measuring the independence of 
subsystems so defined. We then construct a theory based on this definition of subsystem 
together with the additivity of quadratic invariants, and compare its predictions with 
the measured vorticity probability density function. 

1.2. Background 

In a remarkable series of papers Burgers (1929a, b,c, 1933a, b,c,d) [reprinted in Nieuw- 
stadt & Steketee (1995)] appears to be the first researcher to apply statistical mechanics 
ideas to the description of fluid turbulence. Many basic ideas used by later researchers 
were introduced first by Burgers in these rarely cited papers. Burgers introduced both 
lattice and Fourier models and showed that such models satisfy Liouville's theorem when 
viscosity is neglected. He used a counting argument to derive an entropy expression and 
obtained a corresponding entropy maximization principle. He proposed a microscopic 
scale for describing turbulent motion during short intervals of time and defined macro- 
scopic quantities by counting possible streamfunction realizations for sequences of time 
intervals. His analysis is based on the Reynolds stress equation, and he obtained a prob- 
ability distribution that can be used to calculate the mean value of the Reynolds stress. 

Motivated by the work of Burgers, Onsager (1949) took up the subject and considered 
a representation of the vorticity field in terms of a set of point vortices, zero-area vortices, 
of equal strength. Because this results in a finite-dimensional particle-like Hamiltonian 
system, Onsager could proceed to apply techniques of classical statistical mechanics. He 
gave arguments for the existence of negative temperatures and the occurrence of coherent 
structures in a confined region, which are often observed in nature. Related ideas have 
been further pursued by many researchers [e.g. Joyce & Montgomery (1973); Matthaeus, 
Stirbling, Martinez, Oughton & Montgomery (1991); Eyink & Spohn (1993); Yin, 
Clercx & Montgomery (2004)] [see Eyink & Sreenivasan (2005) for a recent review]. For 
example, Joyce & Montgomery (1973) studied the statistical mechanics of point vortices 
within a mean field approximation, and argued that in the negative temperature regime, 
large like-signed vortices are the most probable state. 

T.D. Lee (1952) projected three-dimensional fluid equations (including MHD) onto 
a Fourier basis and truncated to obtain a finite-dimensional system. Evidently unaware 
of the early work of Burgers (1933d), he again demonstrated that his truncated system 
satisfies a version of Liouville's theorem and was thus amenable to techniques of statis- 
tical mechanics. Later, Kraichnan considered two-dimensional fluids [Kraichnan (1967, 
1975); Kraichnan & Montgomery (1980)] and noted that out of the infinite number of 
invariants, two quadratic invariants, the so-called rugged invariants, remained invariants 
after truncation. Kraichnan and Montgomery argued that these rugged invariants are the 
important ones, and obtained an equilibrium state, which is related to that obtained by 
minimum enstrophy arguments put forth by selective decay hypotheses [Leith (1984); 
Maxworthy (1984); Bouchet & Sommeria (2002)]. Also, using Kolmogorov-like dimen- 
sional arguments and the rugged invariants, Kraichnan argued for the existence of direct 
and inverse cascades for two-dimensional turbulence [Kraichnan (1967)]. 

The two-dimensional Euler equation, like the Vlasov and other transport equations, can 
be viewed as mean field theory. Such equations are known to generate fine structure in the 
course of evolution. This led Lynden-Bell (1967) to consider a coarse graining procedure 
coupled with the idea of preserving all of the infinity of invariants such theories possess. 
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He applied his ideas in the context of stellar dynamics, but the ideas are akin to those 
used in treatments of the classical electron gas by generalizations of Debye-Hiickle theory 
[e.g. Van Kampen & Felderhof (1967)]. Later, such ideas were used in the fluid context by 
Robert (1991), Robert & Sommcria (1991), Robert & Sommeria (1992), Miller (1990), 
and Miller, Weichman & Cross (1992), and again in the stellar dynamics context by 
Chavanis, Sommeria & Robert (1996). Our development to a large degree parallels that 
of these authors. In these works a microscopic probability distribution represents a local 
description of the small-scale fluctuations of microscopic vortices. The streamfunction is 
assumed to be uniform on the microscopic scale, and an equilibrium state is obtained 
by maximizing the Boltzmann entropy of microstates, an entropy that is obtained by a 
counting argument first given by Lyndcn-Bcll. This produces a most probable state. 

More recently, the necessity of incorporating the infinite number of invariants in sta- 
tistical mechanics theories has been brought into question, and theories based on finite- 
dimensional models with fewer constraints have been developed. Majda & Holen (1997) 
have argued that including an infinite number of invariants provides no additional sta- 
tistical information, and Turkington (1999) has argued that previous theories have not 
properly handled the neglected small scale phenomena, and he has proposed a theory that 
uses inequality constraints associated with only the convex invariants. Our approach is 
perhaps most closely aligned to these works, but is distinguished by the fact that the in- 
variants chosen are explicitly based on the additivity argument, the choice of subsystems, 
and experimental observation. 

Natural phenomena in atmospheres and oceans have served as a motivation for the 
application of statistical mechanics to two-dimensional fluid flow [e.g. Salmon, Holloway 
& Henderschott (1976)]. Examples include zonal flows in planets, such as the jet stream 
and the polar night jet, and organized coherent vortices, such as the Great Red Spot of 
Jupiter [Maxworthy (1984); Sommeria, Meyers, & Swinney (1991, 1988); Marcus (1993); 
Bouchet & Sommeria (2002)]. Attempts have been made to explain such naturally occur- 
ing phenomena in terms of the coherent structures found to emerge in quasi-geostrophic 
and two-dimensional turbulence after long time evolution. With external small-scale forc- 
ing a few long-lived and large structures resulting from nonlinear merging processes are 
seen to be stable self-organized states that persist in a strongly turbulent environment 
[Mc Williams (1984); Boucher, Ellis & Turkington (2000)]. These structures have been 
studied over many years, often because of their relevance to large-scale geophysical and 
astrophysical flows [Marcus (1993)]. In statistical mechanics, such steady states with 
large structures are envisioned to be the most probable state arising from some extrem- 
ization principle. Various cxtremization principles [e.g. Leith (1984)] have been proposed 
with selected global invariants of the system used as constraints. Observations of tur- 
bulent flow with large coherent structures in a rotating annulus [Sommeria, Meyers, & 
Swinney (1988); Baroud, Plapp, She & Swinney (2002); Baroud, Plapp, Swinney & She 
(2003); Aubert, Jung & Swinney (2002); Jung, Storey, Aubert & Swinney (2004)] have 
led us to reconsider statistical mechanics in the context of rapidly rotating systems. 

1.3. Notation & Organization 
By necessity this paper contains much notation. To aid the reader we give a brief summary 
here. As noted above, statistical mechanics deals with two scales: the microscopic scale 
S, characteristic of microscopic m-cells, and the macroscopic scale A, characteristic of 
macroscopic M-cells. Several averages are considered. The symbol (-)s denotes an average 
with probability density P$, where choices for the subscript S will be used to delineate 
between different cases. The appropriate volume measure will be clear from context but 
is also revealed by the argument of P$. Averages with uniform density are denoted 
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by -< ■ >-s, where the subscript denotes the integration variable. An exception is the 
time average, which we denote by an overbar. Thus, the time average of a function 
is denoted by /, and f = J fdt/T =-< / >~ t . The limits of integration for this kind 
of average will either be stated or will be clear from context. We denote the potential 
vorticity field by q{x, y, t), by which we always mean a function. For the potential vorticity 
distribution on a M-cell (subsystem) we use (, an independent variable. Another source 
of possible confusion is that the symbol (3 is used for the energy Lagrange multiplier, as is 
conventional in statistical mechanics, while the beta-effect of geophysical fluid dynamics 
is embodied here in the symbol h. 

The paper is organized as follows. The experiment is described in §2 and equations that 
govern the dominant physics are reviewed in §3. In §4 we describe some basic ideas about 
statistical mechanics, as needed for the application to the fluid system of interest. In §5 we 
describe statistical mechanics in the mean field approximation and compare predictions 
with experiments. Here we show that predictions of the theory are in accordance with 
experiments. Finally, in §6 we conclude. 



2. Experiment 

The experiments are conducted in a rotating annulus (Fig. 1). The annulus has an 
inner radius = 10.8 cm, outer radius r Q = 43.2 cm, a sloping bottom, and a flat 
transparent lid. The bottom depth varies from 17.1 cm at the inner radius to 20.3 cm 
at the outer radius, giving a bottom slope of n = —0.1. For the data analyzed in this 
paper, the rotation frequency of the annulus is tt/2ir = 1.75 Hz. An azimuthal jet is 
generated in the annulus by pumping water in a closed circuit through two concentric 
rings of holes at the bottom. Fluid is pumped into the annulus through an inner ring at 
r = 18.9 cm and extracted through an outer ring at r — 35.1 cm; both rings have 120 
circular holes. Each hole has a diameter of 2.5 mm, and the total pumping rate is 150 
cm 3 /s. The action of the Coriolis force on the outward flux generates a counter- rotating 
azimuthal jet. A counter-rotating flow is generally more unstable than a co-rotating flow 
[Sommeria, Meyers, & Swinney (1991)]. 

The water is seeded with neutrally buoyant particles (polystyrene spheres, diameter 
150 — 200 /im). Light emitting diodes produce a 3 cm thick horizontal sheet of light that 
illuminates the annulus at mid-depth. The particles suspended in the water are imaged 
with a camera located 2 m above the annulus, and the camera rotates with the tank. 
Particle Image Velocimctry (PIV) is used to obtain the full two-dimensional velocity field 
[Baroud, Plapp, Swinney & She (2003)]. 

The flow can be characterized by the Reynolds, Rossby, and Ekman numbers. The 
maximum velocity £7 max ~ 22 cm/s, the length L = 16.2 cm (taken to be the dis- 
tance between the two forcing rings) and the kinematic viscosity v = 0.01 cm 2 /s yield a 
Reynolds number UL/v = 3.5 x 10 4 , indicating that the flow is turbulent. The Rossby 
number (ratio of incrtial to Coriolis force) is uj rms /2£l = 0.11 (where uj rms is the rms 
vorticity), which indicates that the Coriolis force is dominant, as is the case for planetary 
flows on large length scales. Finally, the small Ekman number, v/2L 2 Q, = 3 x 10~ 4 , indi- 
cates that dissipation in the bulk is small. The Ekman time, te = Li l /2(vfl) 1 ^ 2 (where 
Lh is the mean fluid height) for dissipation in the boundary layers is 30 sec, a time much 
longer than the typical vortex turnover time, 2 sec. The dimensionless numbers indicate 
that the flow is quasi-geostrophic; previous studies of turbulence in the annulus have 
indeed confirmed the strong two-dimensionality of the flow [Baroud, Plapp, Swinney & 
She (2003)]. 
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radius (cm) radius (cm) 

Figure 1. (a) Schematic diagram of the experimental apparatus. The tank rotates at 1.75 Hz. 
Flow is produced by pumping water through a ring of inlets (I) and outlets (0) in the bottom of 
the tank. The Coriolis force acts on the radially pumped fluid to produce a counter-rotating jet. 
(b) The vorticity field and contours of streamfunction at mid-height of the tank, determined from 
Particle Image Velocimetry measurements. The streamfunction contours are equally spaced in 
streamfunction value, (c) The azimuthal velocity averaged over both time and azimuthal angle, 
as a function of radial position, (d) The vorticity (solid line) and streamfunction (dashed line) 
averaged over time and azimuthal angle, as a function of radial position. 



3. Dynamics 

The barotropic assumption is widely used to describe the flow inside the tank. The 
equation of motion for a barotropic fluid with topography is given by 

^ + u-Vq = D + F, (3.1) 

where q = (— V 2 V> + 2fi)/£/, is the potential vorticity, Lh is the tank depth, ip is the 
streamfunction, u = (dip/dy, —dip/dx), D denotes dissipation, such as that due to molec- 
ular viscosity, i/V 2 w, or Ekman drag, —lu/te, and F denotes a vorticity source due to 
the pumping. Often the potential vorticity is approximated by 

q = -W 2 1p + h, (3.2) 

where h accounts for the beta-effect and is here a linear function of radius, h = l^lnrjhh 
where r\ is the bottom slope. Over the years strong evidence has accumulted that (3.1) 
describes the dominant features of the experiment [Sommeria, Meyers, & Swinney (1991); 



Statistical mechanics of two-dimensional turbulence 



7 



del-Castillo-Negrete & Morrison (1992); Meyers, Sommeria & Swinney (1993); Solomon, 
Holloway, & Swinney (1993)]. 

For inviscid flow with zero Rossby number, there is no vertical variation in the velocity 
[Rossby (1939)], and there is evidence that to leading order the drag and forcing terms 
cancel. We are primarily interested in the statistics of motions that occur on the vortex 
turnover time, and these are governed by the inviscid equation, 

^+u-V<z = 0, (3.3) 

which is a Hamiltonian theory. 

A manifestation of the Hamiltonian nature of two-dimensional Euler-like flows such 
as (3.3) is the finite-dimensional Hamiltonian description of point vortices provided by 
Kirchhoff (1883), which played an essential motivating role in Onsager's theory [e.g. 
Eyink & Sreenivasan (2005)]. For a distributed vorticity variable such as q the Hamilto- 
nian form is infinite-dimensional and is given in terms of a noncanonical Poisson bracket 
as follows: 

f t ={q,H} = [i,, q ], (3.4) 

where the Hamiltonian H[q}= j ip(q — h)dxdy/2, and the noncanonical Poisson bracket 
is given by 

~5F SG] 



{F,G} 



dxdy , (3.5) 



Sq ' Sq 

with F and G being functionals, SF/Sq the functional derivative, and [/, g] — f x gy — fy9x- 
Observe that u • Vq — — [ip, q\. This Hamiltonian formulation of the two-dimensional 
Eulcr equation appeared in Morrison (1981, 1982), based on the identical structure for 
the Vlasov-Poisson system [Morrison (1980)], and in Olver (1982). A review of this 
and other formulations can be found in Morrison (1998). The infinite family of Casimir 
invariants C[q] = J C(q)dxdy, where C is arbitrary, satisfies {F, C} — for all functionals 
F, and is thus conserved by (3.3). The presence of these invariants is one way that the 
statistical mechanics of fluids differs from that of particle systems. 



4. Statistical Mechanics & Fluid Mechanics 

As noted in §1, many attempts have been made to apply statistical mechanics to fluids 
and other infinite-dimensional systems. In this section we introduce our notation and 
discuss some basic ideas. 

4.1. State variables 

In classical statistical mechanics the microscopic dynamics is governed by Hamilton's 
equations and the phase space is the 2iV-dimensional manifold with canonical coordi- 
nates (Q a , P a ), Oi = 1, 2, . . . , N, where (Q\, . . . Qn) is the configuration coordinate and 
(Pi, . . . , Pn) is the corresponding canonical momentum. Typically N, the number of de- 
grees of freedom, is a very large number ~ 10 23 . We call this 2A^-dimensional phase space 
T, a standard notation introduced by P. and T. Ehrenfest [Ehrenfest & Ehrenfest (1959)]. 
Our fluid is assumed to be governed by (3.3), an infinite-dimensional Hamiltonian the- 
ory, and thus the instantaneous state of our system is determined by the vorticity-like 
variable q(x,y), which we suppose is contained in some space of functions Q. The index 
a for coordinates of V is analogous to the Eulerian position (x,y), a point in the physical 
domain occupied by the fluid, which is viewed as an index for Q. 

In conventional statistical mechanics, the microscopic dynamics is finite dimensional, 
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and one attempts to explain phenomena on the macroscopic level by considering the ther- 
modynamic limit in which N — > oo. However, for a fluid, the dynamics is already infinite 
dimensional, and thus as noted in §1, to apply statistical mechanics researchers have intro- 
duced various finite-dimensional discretizations. Onsager's description of the continuum 
vortex dynamics in terms of a collection of point vortices amounts to the specification of 
the coordinates of the manifold analogous to T as the spatial positions of the point vor- 
tices, {x\, . . . xn, 2/1, ... , y^). Alternatively, Lee's representation of a three-dimensional 
fluid in terms of a truncated Fourier series has the Fourier amplitudes being coordinates 
of a space analogous to Y. This procedure was carried over to two dimensions by Kraich- 
nan and Montgomery [Kraichnan & Montgomery (1980)]. For our potential vorticity 
variable the Fourier amplitudes are given by = J expi(k x x + k y y) q(x, y)dxdy, where 
k = [k x , k y ). Another alternative is to replace the continuum vorticity by a lattice model 
[e.g. Burgers (1929a); Robert (1991); Robert & Sommeria (1991, 1992); Miller (1990); 
Miller, Weichman & Cross (1992); Majda & Holcn (1997); Turkington (1999)], i.e., an 
expansion in terms of tent functions or finite elements of scale size 5. In the present con- 
text the vorticity is replaced by its values on the lattice, q% = J Ki(x, y; x il y i )q{x 1 y)dxdy, 
where the kernel K t is typically chosen to represent a square lattice with a finite number 
N of sites located at (x i7 yi). In general N — N x N y , where N x and N y are the number of 
lattice points in the x and y directions, respectively. We will refer to this discretization 
as a division into m-cells. 

Given a finite-dimensional system one can make various assumptions, e.g., the proba- 
bilistic assumptions of 'molecular chaos', but this requires a notion of phase space volume 
conservation. 

4.2. Phase Space Volume & Liouville's Theorem 

In classical statistical mechanics one calculates averages over the manifold T, and the 
natural volume element is given by U N a=1 dQ a dP a . However, for Q the situation is not 
so straightforward, and so we explore candidates for the analogous volume element. 

4.2.1. Volume element 

The calculation of averages in a statistical theory requires a phase space measure, 
7S>q, which is a sort of volume element for Q. The volume element can be interpreted 
as a probability measure defined on functions that take values between q and q + dq. 
Averages calculated using the probability measure are functional integrals akin to those 
used in Feynman's path integral formulation of quantum mechanics and in field theory 
[e.g. Schulman (1981); Sundermeyer (1982)]. The various discretizations introduced 
above have been employed to give meaning to functional integrals, but the Fourier and 
lattice models are most common. 

For the Fourier discretization, Kraichnan and Montgomery used the volume element 
T>q = Jlk^k' where the product is truncated at some maximum wave number. Alter- 
natively, the volume element for lattice models is written as Vq = dqi , where dqi is 
a volume element associated with the potential vorticity varying from q to q + dq in a 
lattice partition (xi,yi), and N = N x N y is, as above, the number of lattice sites, which 
have a scale 8. Here, a total volume element Vq is a product of volume elements of each 
lattice site dqi. In the case of a finite small lattice, dqi becomes a one-dimensional vol- 
ume, i.e., dqi = q(xi, yi) + dq(xt, y{) — q{xi, y{) at the lattice point (xi, yi) of the physical 
two-dimensional space. In order for a notion of measure based on phase space volume to 
be useful, the volume must be preserved in the course of time. 
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4.2.2. Liouville's theorem 

Preservation of phase space volume is assured by Liouville's theorem, an important 
theorem of mechanics. As noted above, Burgers and Lee showed that a version of Liou- 
ville's theorem applies to the system governing the Fourier amplitudes for the inviscid 
fluid. For vorticity dynamics the amplitudes satisfy 

^ k = J2^W~( qi ~ fe i)9m, (4.1) 

l.m 

where h\ is the Fourier transformation of the beta effect and £kim = z-(lx m)<5(k+l+m) is 
completely antisymmetric, i.e., ekim = — eikm = — £mik and ekkm = £kik = 0. Therefore, 
antisymmetry directly implies Liouville's theorem, ^^dq^/dq^ = 0. 

Similarly, we have shown directly that the lattice model possesses a version of Liou- 
ville's theorem, which we recently discovered was anticipated in Burgers (1929b). This 
result was also inferred in Turkington (1999). We assume periodic boundary conditions. 
The lattice model discretization can be viewed as an expansion of the vorticity in terms 
of a tent function basis [e.g. Fletcher (1980)]. Upon multiplying the equation of motion 
by a basis function, as is typical with Galerkin projection, and representing all derivatives 
as differences, Eq. (3.3) becomes 

Qi = y~] B ijk ipjq k , (4.2) 

j,k 

which is an equation for the potential vorticity at the lattice point i. Assuming a periodic 
lattice, the quantity Bij k is easily seen to be completely antisymmetric, i.e., B^ k — 
—Bjik = —Bkji and B^i = Buk = 0, just as was the case for ekim- Therefore, Liouville's 
theorem follows, 

E || = E B vk(M jiqk + S k ^) = Yjn'„ Al + Bm^) = , (4.3) 

where each term of the last sum vanishes. Here the matrix M represents the inverse 
Laplacian and B' is another matrix that has the same antisymmetry property as B. 

4.3. Canonical equilibrium distribution 

Having defined phase space and verified Liouville's theorem, we are poised to write a 
partition function and to define phase space averages. The natural expression for the 
partition function associated with the canonical (Gibbs) ensemble is 

Z c = [ e-^W-eWVq, (4.4) 
Jg 

where H is the Hamiltonian of §3 and C denotes the infinite family of Casimir invariants. 
Averages corresponding to (4.4) are given by 

{F) c = f F[q}P c [q;p,C}Vq, (4.5) 
Jg 

where F is a functional of q and the phase space probability density is given by 

P c [q;(3,C] = Z^e^M-CM _ (4 6) 

Expressions (4.4) and (4.5) are functional integrals [Schulman (1981); Sundcrmeyer 
(1982)], and the intent is to give them meaning by discretizing as in §§4.1 and then 
taking the limit N — > oo and 5 — > 0. Finding unique well-defined results with this 
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procedure for such integrals, with other than quadratic functionals in the exponent, is 
usually a difficult task. Consequently, a mean field approach has been taken, which we 
turn to in §5. 

An alternative to the direct evalution of (4.5) is to appeal to the fact that the dy- 
namics of (3.3) is an area-preserving rearrangement [e.g. Lieb & Loss (2001)]. This 
means for an initial condition q , the solution at time t is given formally by q(x,y,t) = 
q (x (x,y,t),y (x,y,t)), where (x (x,y,t),y (x,y,t)) are the initial conditions of the 
characteristics, which satisfy d(x ,y )/d(x,y) = 1. The Casimir invariants are associ- 
ated with relabelling symmetry [e.g. Salmon (1982); Padhyc & Morrison (1996)] and 
possess the same value when evaluated on functions that are related by rearrangement. 
Thus, if one restricts the domain of integration Q to be rearrangements of a given func- 
tion, denoted by Qn, then we should obtain the same answer because {F[q])n = F[q] for 
functionals with integrands that depend only on q, such as Casimirs and exp(C[q]). Here 
( )tz is defined with Pn[q] P] = -Z^ 1 exp (— (3H[q\) and Z-r = / exp (—f3H[q])T>q. 



5. Mean Field Approximation & Statistical Independence 

It is well-known that vorticity equations like (3.3), the Vlasov equation, and other 
transport equations develop fine structure in the course of time. Because of this Lynden- 
Bell (1967) proposed a coarse graining procedure to obtain a most probable state. He 
divided phase space up into hyper-fine cells that are assumed to be capable of resolving 
the fine structure. These are the m-cells referred to in §§4.1, which have a scale size S. 
Experimentally 5 is determined by the resolution, but in ideal theory the fine structure 
can become arbitrarily fine and so a limiting procedure is required. In addition Lynden- 
Bell (1967) proposed larger cells, which we have called M-cells, that characterize a 
macroscopic scale A. The M-cells contain many m-cells that can be freely exchanged 
within an M-cell without changing any macroscopic quantity. Thus one is able to count 
states and obtain an expression for a coarse grained or mean field entropy that can be 
maximized subject to constraints. Later, Miller (1990) and Robert & Sommeria (1991) 
reconstructed and improved this formulation. Miller defined m-cells and M-cells based 
on scales with the property that the energy averaged over M-cells approximates the 
energy averaged over m-cells. However, we argue that the most important condition for 
separating the M-cell and m-cell scale lengths is statistical independence, which assures 
near independence of the probability densities of M-cells, which are viewed as subsystems, 
and is associated with near additivity of the constraints. These are crucial properties. 

Experimentally the two scales can be demonstrated as in Fig. 2. Observe in the upper 
plot of this figure the fine scale structure in the potential vorticity, while in the lower 
plot the streamfunction, due to the integration over the Green's function, is considerably 
smoother. We take the upper scale to be 5 and the lower scale to be A. 

5.1. Counting states 

According to Lynden-Bell's statistics, the number of ways to distribute m-cells into M- 
cells is 

jV r ' -i-r N^l 

w = \ [ n^. [ t(NU-Xrrt I) y (5 ' 1} 

where N r is the total number of m-cells with the rth value of potential vorticity in the 
whole space, and N^ 1 ^ is the total number of m-cells with the rth value of potential 
vorticity in the 7th M-cell. Also, is the total number of m-cells in the Jth M- 

cell. The first product in Eq. (5.1) represents the number of ways to distribute N r 




tt 3tt/2 2tt 

e 



Figure 2. The time-averaged potential vorticity (top two figures) and the streamfunction (bot- 
tom two figures) in the Rossby wave frame. The figures on the left show the fields in the ro- 
tating tank; the figures on the right show the same fields unwrapped. The streamfunction field 
is smoother than the potential vorticity field since the vorticity is given by a second derivative 
(the Laplacian) of the streamfunction. Hence the characteristic length scale in the azimuthal 
direction is larger for the streamfunction than for the potential vorticity. 



m-cells into groups of {N^}, where / counts all M-cells and the second product is the 
number of ways to distribute inside an M-cell. Also, — J2 r -Wr can be understood 
as the number of empty m-cells. Lynden-Bell proposed this manner of counting for stellar 
dynamics [Lynden-Bell (1967), see also Chavanis, Sommeria & Robert (1996)], where 
m-cells represent stars, which are considered to be distinguishable, and there may be 
empty m-cells. However, the statistics for the two-dimensional continuum Euler model 
is a special case of Lynden-Bell's general counting procedure. 

In Miller's application to the two-dimensional continuum Euler model, he assumes that 
all m-cells are occupied by a vortex and these vortices are indistinguishable if they have 
the same value of vorticity. Because there are no empty m-cells, — J2 r and 

because the m-cells are indistinguishable, a factor of J\ r l/(N r l) is added. This counting 
produces 

W = IT T7T • ( 5 - 2 ) 

The above equation already involves statistical independence among different M-cells. 

Boltzmann articulated the entropy as a measure of the number of possible configura- 
tions of the system. Therefore, the entropy S is defined to be the logarithm of the total 
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number of configurations, InW. If is large, Stirling's formula gives 

5 = W,-^(^))ln^). (5.3) 

In the continuum limit of potential vorticity levels, /N^ is replaced by Pm((', x, y), 
and J2 r i I d(dxdy. In short, the index / represents the coordinates for the discretized 
M-cells and the index r represents the ordered level sets of potential vorticity inside the 
M-cells. Thus, it is replaced by the continuum vorticity variable (, the vorticity on an 
M-cell. With these observations, the resulting total mean field entropy is seen to be 

[Pm] = - J Pm((; x, y) In P M ((;x,y)d(dxdy = - J (In P m )m dxdy (5.4) 

where Pm{(] x, y) is the probability density in the mean field approximation. The density 
of Pm(Ct x iU) is centered at the point (x,y) and satisfies the normalization /Pm^C = 1- 
The integration over dxdy can be viewed as a sum over the M-cells that cover the domain 
of the fluid. The second equality of (5.4) follows from the definition (A)m = J APMd(, 
and thus Sm[Pm] can be naturally termed the (mean field) Boltzmann-Gibbs entropy. 

In closing this subsection, we reiterate that the potential vorticity variable q is a 
field variable, a function of coordinates. However, when we introduced the probability 
density Pm on M-cells, we used £, an independent variable, to represent the values of 
the potential vorticity on an M-cell. 

5.2. Mean field canonical distribution 
Given the mean field entropy Sm we can proceed to obtain the mean field density 
Pm{(\ x, y) as the most probable state by extremization subject to particular mean field 
constraints. These constraints and their corresponding Langrange multipliers are given 
as follows: 

(a) The Hamiltonian constraint is obtained by replacing the vorticity variable q in 
H[q] with its mean field average, to obtain a mean field energy, 

H m [Pm] = \f ((P M (C,x,y)-h) CPM((';x',y')G(x,y;x',y')d(dxdyd{'dx'dy' 

m ((C) m - h) dxdy (5.5) 



\ J M> 



where (()m = J (Pm d( and (V')m is defined by 

(C)M = -V 2 W M + fc. (5.6) 

The Lagrange multiplier associated with this constraint is taken to be the constant value, 
— (3, where the minus sign is by convention. 

(b) The normalization constraint is / Pm^C = 1- This is a normalization on each M- 
cell; thus, although Pm depends on position, the integration does not. Because this is a 
constraint for each point (x, y), the Lagrange multiplier in this case depends on position. 
We call it j(x,y), and the quantity that appears in the variational principle is 



N 



[P M } = J \(x,y)P M {<;;x,y)d(dxdy. (5.7) 

(c) The mean field Casimir constraint, roughly speaking, contains the information that 
on average, the area between any two contours of vorticity remains constant in time. More 
precisely, the quantity g(() = J Pm dxdy is taken to be constant. Because this is true 
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for all C) the Lagrange multipler n is likewise a function of £ and the constraint can be 
written as 

C m [Pm] = -fij KCMOdC = -/3 J »(OP M ((;x,y)d(dxdy, (5.8) 

where the prefactor of —0 is again by convention. This constraint is the mean field version 
of the family of Casimir invariants C[q]. 

Now we are in position to obtain the most probable state by extremizing the quantity 
Fm = Sm — (3Hm + Nm + Cm, i.e., upon functional differentiation with respect to Pm, 
SFm/SPm — implies 

P M (C,x,y;f3,i,) = Z^e-MW™-^ , (5.9) 

where Zm — J e~^^ M ~ ,l ^dC, and evidently Pm is normalized. Equation (5.9) is 
the mean field counterpart to (4.6) and could aptly be termed the canonical (Gibbs) 
mean field distribution. The above variational principle and extremal distribution (5.9) 
appeared in essence in an appendix of Lynden-Bell (1967). 

Given (5.9) we are in a position to calculate (C)m and then substitute the result into 
(5.6). This gives the mean field Poisson equation, 

V 2 (V>)m = Z M l j Qe-^WM-n(i)) d Q + h . (510 ) 

Versions of this equation have been solved in various references [e.g. Robert & Sommeria 
(1991); Miller (1990); Majda & Holen (1997)], but we will not do this here. 

We conclude this subsection by giving a heuristic connection between ( )m, a prescrip- 
tion for averaging functions, and ( ) c , a prescription for averaging functionals. Consider 
the functional q[x! , y'), by which we mean the evaluation of the function q at the point 
(x',y'), and evaluate 

(q(x',y')) c = [ q(x',y')P c [q;p,C}Vq. (5.11) 
Jg 

If we rewrite (5.11) as an integral on M-cells, where q(x', y') is qp, write Vq = J\j dqj, 
and then assume statistical independence of M-cells, P c = flj Pi, we obtain 

(q(x',y')) c = J qp Y[Pr \\dq,, = J q v P P dq P = J (P M d( = (Q M . (5.12) 



This derivation emphasizes the need for near statistical independence of M-cell subsys- 
tems. 

5.3. Ruggedness & additivity 

Classical statistical mechanical treatments of the canonical ensemble allow for subsys- 
tems to interact and exchange energy, but their interaction is assumed to be weak and the 
details of the interaction are usually ignored in calculations. Neglect of the interaction 
energy results in the energy being equal to the sum of the energies of the individual sub- 
systems, i.e., the energy is an additive quantity. In conventional treatments only additive 
invariants are used in calculating the most probable distributions, and in some treat- 
ments [e.g. Landau & Lifshitz (1980)] this requirement is explicitly stated. The reason 
for this is that additive invariants give rise to statistical independence of subsystems. In 
our treatment of fluids, subsystems are M-cells and so we consider invariants that are 
additive over these regions. There is a close connection between ruggedness of invariants 
and the property of additivty. We show that only the rugged invariants are additive, and 
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thus they characterize the statistical properties of M-cells. In §§5.4 and §§5.5 we will see 
that experimental results support this reasoning. 

Kraichnan & Montgomery (1980) Fourier transformed and truncated to obtain a 
finite-dimensional system. They argued that the truncated remnants of the total vortic- 
ity, enstrophy, and energy are the only invariants to be used in a statistical mechanics 
treatment because these invariants are rugged, i.e., they remain invariants of the trun- 
cated system. They also appear to be aware that these invariants possess the property of 
additivity, but they do not emphasize this point. Although Turkington (1999) has argued 
that this kind of truncation does not properly handle small scale behavior, we find that 
this theory does a fairly good job at predicting the energy spectrum, but we will report 
on this elsewhere. We argue in general that such invariants are important because they 
are the only additive invariants. Below we consider a somewhat more general setting. 

Because of Parseval's identity, the quadratic invariants are additive and higher order 
invariants are not. To see this, suppose we define M-cells to be composed of amplitudes 
of some subsets of Fourier modes, which we denote by m. Then a sum over modes can 
be done in groupings, i.e. ^ k = J2i J2 KI ■ (This is the idea behind spectral reduction 
[Bowman, Shadwick & Morrison (1999)], a computational method where groupings of 
Fourier modes (bins) are described by a single representative.) For the quadratic Casimir 
invariant, the enstrophy, we have 



Co 



f q 2 dxdy=(2n) 2 Y / \^\ 2 , (5-13) 



and defining an M-cell enstrophy by = (27r) 2 J2 Ki \qk\ 2 , we obtain C 2 = ^ZjC^- 
Similarly, the energy can be written as a sum over M-cell energies, 
linear Casimir invariant C\ = J qdxdy merely reduces to the zeroth Fourier coefficient, 
and is thus in a trivial sense additive. Higher order invariants, C n = J q n dxdy for n > 2, 
have Fourier representations that are not reducible to expressions in terms of a single 
sum over M-cells. 

The discretized lattice model has properties similar to those described above. The 
quadratic Casimir invariant and energy reduce to sums over a finite number of m-cell 
lattice variables, qi, hi and ipi, which are potential vorticity, height, and streamfunction 
represented in terms of the kernel function Ki of §§4 as follows: 

C 2 = J q 2 dxdy = E / K i K j ( H ( lo dxdy = E ''■ / >i'li > 
H = \ J qndxdy = \ ^2 f K * K k(l* - h i)^k dxdy 

i.k 

= ^ y^(g» - hi)Z ik i[)k = y^ili ~ h i)ZiMi - h i) ( 5 T4) 

i.k i.j 

where Zij = J KiKj dxdy and Zij = ^2 k ZikMkj are symmetric commuting matrices. 
These invariants are rugged, i.e, they are conserved by the finite dynamical system ob- 
tained by projection onto the lattice. In addition, because Z and Z commute, one can 
always find an orthogonal matrix O that satisfies Z = O^O and Z = O^O, where 
Dij = diSij and Dij = diSij are diagonal matrices. Defining q 1 = qO, h! = hO and 
ip' = ipO, the enstrophy and energy become 
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(5.15) 




where 7 is the index for the 7th M-cell and Ki denotes the set of m-cells in the 7th 
M-cell. 

This coordinate transformation simultaneously diagonalizes the quadratic Casimir in- 
variant and the energy. However, higher-order Casimir invariants are in general not 
rugged and are in general not simultaneously diagonalizable. Thus, higher order invari- 
ants are not additive, which means M-cells share contributions from these invariants. 
In this sense, invariants of order higher than quadratic arc not useful for describing the 
statistics of M-cells, which by assumption are independent. 



Now we turn to the question of how to find subsystems, i.e, how to a find a good defini- 
tion of the M-cells. First we note that flows inside the rotating tank with the sloped bot- 
tom have azimuthal undulations in most physical quantities (streamfunction, potential 
vorticity, etc.), and these undulations have been identified as Rossby waves [del-Castillo- 
Negrete & Morrison (1992); Solomon, Holloway, & Swinney (1993)]. In a co-rotating 
frame, these waves propagate in the rotation direction at constant velocity. Thus, by 
shifting to a frame moving at the phase velocity of the Rossby wave, we obtain a pattern 
that is statistically stationary on large scales. For example, the wavy patterns correspond- 
ing to the time-averaged streamfunction and potential vorticity are shown in Fig. 2. As 
noted before, the streamfunction is fairly smooth, characteristic of the scale A, is mono- 
tonically decreasing in the radial direction, and describes a strong zonal flow. However, 
the time-averaged potential vorticity is scattered with fine structure in space, the 5 scale, 
but still has a wavy mean pattern similar to that of the time-averaged streamfunction. 
So, this suggests that the first step toward defining M-cells is to consider a frame moving 
at the phase velocity of the Rossby wave. 

Having determined the frame, we seek M-cells that are statistically independent. Be- 
cause strong correlation in a preferential direction might affect the geometry of M-cells 
and associated additive invariants, we have measured the correlation function, 



where (9, r) are the usual polar coordinates. From a large data set of PIV measurements 
we obtain the time average of the velocity field, whence we calculate the potential vorticity 
at different positions. Then the integrals of (5.16) are performed with the spatial limits 
being the bulk of the area occupied by the fluid with a resolution of S « 0.8 cm and the 
time limit taken to be 80 revolutions with 47 measurements. The result of this procedure 
is presented in Fig. 3, which shows contours of C cor plotted on a AO— Ar plane. The highly 
anisotropic nature of the contours suggests there is significantly less correlation in the 
radial direction than in the azimuthal direction. Thus to achieve consistent independence 
the shape of an M-cell should be elongated. 

In the course of tracking blobs of fluid we generally observe that to good approximation 
such blobs follow contours of the time-averaged streamfunction. This, together with the 
the C cor plot, suggests that a good coordinate for dividing the system into subsystems is 
the time-averaged streamfunction, 



5.4. Statistically independent subsystems 
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Figure 3. Contours of the correlation function C cor (Ar, rAO) illustrating the anisotropic nature 
of the potential vorticity field, which has longer range correlation in the azimuathal direction 
than in the radial direction (cf. Fig. 2) . 



Contours of ?/> tend to be smooth and, we argue, are part of a natural coordinate system 
for describing turbulence with a mean flow that has slow spatial dependence. (We have 
also considered q but found it to be not as good because of its greater variability.) To 
complete the coordinate system, we introduce a coordinate x, which is conjugate to V> 
and therefore satisfies 

1 d^d\ _ 1 dj)_dx = 1 (51g , 
r 09 dr r Or 30 

Thus the coordinate transformation (0, r) < — ► (x, 4>) satisfies rdrdO — dxdO. 

We propose that contours of tp define M-cells, which we take to be of small (infinitesi- 
mal) width in this coordinate, and we propose that the \ coordinate at fixed -tp represents 
a continuum of m-cells. We imagine an M-cell to be a region (nearly a curve) at fixed ip. 
Hence with this definition, the probability density Pm sxp i depends only on the potential 
vorticity variable ( and on the coordinate r/i; i.e., Pm cxv (C> VO is the probability of finding 
a potential vorticity value ( in the tp M-cell. Thus the ensemble average of an arbitrary 
function / is written as 

(f)M ex M= I f(^)PM exp (CA)dC, (5.19) 

where Pm exp is normalized as J PM exp dC, = 1. In practice we can determine the probability 
Pm^ p from data by the relative frequency definition (cf. §§5.5), and then proceed to cal- 
culate (5.19). However, this is equivalent to averaging over x an d t; e.g. (Qm exp = -<9^ x , 
where -< q >- x = / qdxj J d\- Given (C)m cxp and using (5.6) to define (V , )m C!CP we simi- 
larly have the equivalence (ip) M exp {i>) = ~<ip('<P,X]t)^x = V'' where the second equality 
follows by definition. The undular streamfunction of Fig. 2 mainly represents Rossby 
waves. These wavy patterns are quite robust and often behave as barriers to mixing. 
In the Rossby wave frame, our data indicate that the instantaneous streamfunction is 
close to the time-averaged streamfunction, i.e., ~<ip(i>,Xit) deviates from ■tp by less 
than 10 percent. The above comments can be viewed as an experimental verification of 
ergodicity. 

In terms of the above notation the energy and enstrophy densities on M-cells can be 
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Figure 4. (a) Enstrophy fluctuations AtC2(V>) [Eq. (5.21)] and energy fluctuations AtH(^) 
as a function of tp. The fluctuations are small, indicating that energy and enstrophy are nearly 
conserved for our choice of subsystem, (b) Total enstrophy variations A^C2(t) [Eq. (5.22)] and (c) 
total energy variations A^i^i) with time; the variation is small, indicating that the quantities 
for our choice of subsystem are almost conserved in time. 



written as 



{H) Mex M = 



I 



(5.20) 



and two quantities that measure spatial and temporal fluctuations of these invariants 
can be compactly written as follows: 



A*C 2 (i) 



1/2 



1/2 



(5.21) 



(5.22) 



with similar expressions for AtH(iP) and A^H(t). Figure 4 depicts these quantities. Panel 
(a) shows temporal fluctuations as a function of the spatial coordinate The middle 
regions of the experiment, where strong zonal flows exist, is describable by statistical 
mechanics. However, near the walls, corresponding to high and low ip values, statistical 
mechanics fails because of large fluctuations. Similarly, in panel (b) the spatial fluctua- 
tions are plotted versus time, and it is observed that these fluctuations are quite small. 
We have measured similar quantities for the cubic and quartic Casimir invariants and 
the fluctuations are two or three times greater. 

An integrated measure of the goodness of our streamfunction based M-cells is displayed 
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| Fluctuation Measure | 




|^A T C 2 ^| 


j square cells j 


| 0.2233 


| 0.6425 | 


| streamfunction cells | 


| 0.0343 


| 0.0627 | 



Table 1 . Comparison of fluctuations for square cells with our streamfunction based cells. Both 
the energy fluctuation measure -< AtH y$ and enstrophy fluctuation measure -< AtC-2 y^ are 
considerably smaller with the streamfunction based cells. These small fluctuations allow the 
division of the system into M-cells, consistent with the statistical independence and additivity 
assumptions of statistical mechanics. 



in Table 1. Here we have integrated -< AtH y^ and -< AtC2 over central values of 
ip and compared them with counterparts derived using square cells. By this measure 
streamfunction based cells are nearly ten times better than square cells. 

Thus, in summary we have strong evidence supporting the use of streamfunction based 
M-cells. The evidence of Fig. 4 and Tabic 1 imply both statistical independence and the 
additive nature of the quadratic invariants of these macro-cells. 

5.5. Prediction for PDFs 

Based on the arguments in the previous section, we consider only two invariants out of 
the infinitely many invariants conserved by the ideal dynamics. Consequently, we obtain 
the following equilibrium distribution: 

P Mexp (C^)=Z-l p e-^-^\ (5.23) 

where Zm cx = J e~ /3 '^ ^ d( depends only on ip. Note, the function h has cancelled out 
in the normalization. This probability density function (PDF) has the form of Gaussian 
that is shifted by (3ip/2j, a position dependent term that can be interpreted as a sort of 
'local wind'. 

In Fig. 5 we compare (5.23) with experimental results. Figures 5(a) and 5(b) show 
that experimental data on a typical M-cell closely agree with the Gaussian distribution 
of (5.23). Each distribution is shifted by its mean value of potential vorticity (C)M exp - 
Figures 5(c) and 5(d) show the total probability P tot (Q, which is the sum of the prob- 
abilities over all the M-cells, i.e., P tot (() =-<PM exp {(;4>) These plots are decidedly 
non-Gaussian. 

The source of the small deviation in the tails of the distribution plotted in Fig. 5(b) can 
be attributed to experimental resolution and sample size. We analyzed this by considering 
the scaling of the skewness and kurtosis with respect to the number of data points N 
and the experimental subsystem width A-tp. For a gaussian distribution the kurtosis is 
three and the skewness is zero. If we fix N (large), our data indicate that the kurtosis 
approaches three as (Aip) a , where the exponent a is less than unity and is approximately 
0.5 for our last (smallest Aip) data points. At fixed Axp (small) we find that the tails 
decrease as we increase N. A similar examination of the skewness reveals randomness 
about the value zero. 

The next question is, what is the most probable value of potential vorticity in each 
M-cell? The probability distribution of Eq. (5.23) gives a relation between the averaged 
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Figure 5. The measured probability distribution of potential vorticity (data points) on a typical 
M-cell is nearly Gaussian (dashed line), in accord with (5.23), as illustrated by these plots on 
(a) linear and (b) logarithmic scales. In contrast, the potential vorticity of the whole system, 
shown in (c) and (d) respectively, departs significantly from a Gaussian. 



vorticity and the streamfunction, 

(C)m cxp - J OPaw (C; $ d( = -e$ , (5.24) 

which follows by elementary integration. Here e = (3/(2j) is the ratio of two Lagrange 
multipliers. Figure 6 shows a linear relation between the ensemble-averaged potential 
vorticity (C)m„ xp and the time-averaged streamfunction ip, as predicted by Eq. (5.24). 

Therefore, our theoretical predictions based on a mean field approximation are in 
good accord with PDFs on M-cells and the averaged values of potential vorticity and 
streamfunction from experiments. Our theory also indicates that equilibrium can be 
locally achieved in M-cells, even though the system as a whole is turbulent and non- 
Gaussian. 



6. Conclusions 

In this paper we have emphasized the relationship between additive invariants and 
statistical independence: probability densities that result from entropy maximization 
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Figure 6. The ensemble-averaged potential vorticity (g)M cxp exhibits a dependence on the 
time-averaged streamfunction tp that is linear except near the walls (at the ends of the range of 
tp). The dots are mean values of (q)M exp , and the vertical lines correspond to standard deviations 
of (q)M exp at a fixed tp. The data fit well the straight line (a least-squares fit), in accord with 
the prediction of Eq. (5.24), where the slope is the ratio of two Lagrange multipliers. 

principles, such as that of §§5.2, will decompose into a product over subsystems if the 
entropy is logarithmic (extensive) and the invariants included as constraints are additive 
over subsystems (M-cells). We have also emphasized that additivity and, consequently, 
independence depend on the definition of subsystem. This idea appears, at least implicitly, 
in conventional statistical mechanics. For example, in the classical calculation of the 
specific heat of a solid, where one considers a solid to be a collection of lattice sites with 
spring-like nearest neighbor interactions, the Hamiltonian achieves the form of a sum over 
simple harmonic oscillators. However, such a diagonal form requires the use of normal 
coordinates, and only then is the partition function equal to a product over those of the 
individual oscillators. Thus the notions of subsystem (here a single oscillator), additivity, 
and statistical independence are intimately related. 

In our application of statistical mechanics to inhomogeneous damped and driven turbu- 
lence, we have discovered experimentally that a good definition of subsystem is provided 
by the temporal mean of the streamfunction. With this definition, the quadratic invari- 
ants (energy and enstrophy) are additive, and the concomitant probability density of 
(5.23) agrees quite well with experimental results for both the distribution of vorticity, 
as depicted in Fig. 5 (a) and (b), and the mean state, as depicted in Fig. 6. 

An alternative interpretation of our results can be obtained by the counting argument 
of §§5.1. Our definition of subsystem amounts to the idea that potential vortices on the 
same contour of time-averaged streamfunction can exchange their positions with little 
change in the energy and enstrophy. However, the relocation of two potential vortices 
that arc on different contours of the streamfunction should result in a large change of 
the invariants. In this sense, the number of possible configurations in phase space can be 
counted, and the maximization of the entropy so obtained gives our result. 

Our discussion of statistical independence and additivity has been heuristic, in the 
spirit of Boltzmann and Gibbs. We suggest that a more rigorous development could use 




Statistical mechanics of two-dimensional turbulence 



21 



the techniques described in other works [e.g. Miller, Weichman & Cross (1992); Majda 
& Holen (1997); Turkington (1999)] adapted to our r/5 coordinate that describes our 
subsystems. For example, one could begin with an appropriate sequence of lattice models 
and obtain a continuum limit. 

Although in this paper we have focused on a geostrophic fluid, our procedure is of 
general utility and is applicable to physical systems governed by a variety of transport 
equations. The unifying formalism is the noncanonical Hamiltonian description of §3, 
which plays the unifying role played by finite-dimensional canonical Hamiltonian sys- 
tems in conventional statistical mechanics. Thus we expect our approach to apply to 
Vlasov-Poisson dynamics, kinetic theories of stellar dynamics, drift-wave plasma mod- 
els, and other single-field models that possess the noncanonical Poisson bracket of (3.5). 
Generalization to multi-field models such as reduced magnctohydrodynamics, stratified 
fluids, and a variety of physics models governed by generalization of the Poisson bracket 
[Thiffcault & Morrison (2000)] of (3.5) provides an avenue for further research. 
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